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X-ray observations of the Large Magellanic Cloud pulsar 
PSR B0540-69 and its PWN 
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ABSTRACT 

PSR B0540-69 is a young pulsar in the Large Magellanic Cloud that has similar 
properties with respect to the Crab Pulsar, and is embedded in a Pulsar Wind Nebula. 
We have analyzed the complete archival RXTE dataset of observations of this source, 
together with new Swift-XRT and INTEGRAL-IBIS data. Accurate lightcurves are 
produced in various energy bands between 2 and 60 keV, showing no significant energy 
variations of the pulse shape. The spectral analysis shows that the pulsed spectrum is 
curved, and is best fitted up to 100 keV by a log-parabolic model: this strengthens the 
similarities with the Crab pulsar, and is discussed in the light of a phenomenologic 
multicomponcnt model. The total emission from this source is studied, the relative 
contributions of the pulsar and the PWN emission are derived, and discussed in the 
context of other INTEGRAL detected pulsar/PWN systems. 

Key words: stars: pulsars: general - stars: pulsars: individual: (PSR B0540-69) - 
X-rays: stars - gamma rays: observations 
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1 INTRODUCTION 

PSR B0540-69 is situated in the Large Magellanic Cloud, 
and is one of the few known extragalactic pulsars. Its pulsed 
emission was discovered by Seward et al. (1984) in the soft 
X-ray band with the Einstein observatory, and thereafter 
also in the optical (Middleditch & Pennypacker, 1985) and 
radio bands (Manchester et al., 1993). 

This pulsar is usually referred as the "Crab Twin" , be- 
cause of the similar period and period derivative with re- 
spect to the well known Crab Pulsar (PSR B0531+21, see 
Table [T] for a comparison). However, their pulse shapes are 
strongly different: for the Crab it consists of two peaks, with 
a phase separation of about 0.4, that are visible at all wave- 
lengths, while PSR B0540-69 has a single, broadly symmet- 
ric pulse about 0.5 wide in phase, with some structure at 
the maximum. 

This pulsar, like the Crab, is also embedded in a Pul- 
sar Wind Nebula (PWN) that shines in the optical to X-ray 
band (Serafimovich et al., 2004). This pulsar/PWN system 
has been observed many times in the X-rays, by ROSAT 



E-mail address: riccardo.campanaOunironial.it 



Table 1. Comparison between various parameters for 
Crab Pulsar (PSR B0531+21) and PSR B0540-69. 
Data are from the ATNF Pulsar Database, available at 
http : //www . atnf . csiro . au/research/pulsar/psrcat/. 

Crab PSR B0540-69 



Period P (ms) 


33.5 


50.5 


Period derivative P (xlO 13 s/s) 


4.2 


4.8 


Braking index n 


2.5 


2.1 


Distance d (kpc) 


2 


50 


Characteristic age r (yr) 


1240 


1670 


Magnetic field B s (xlO 12 G) 


3.78 


4.96 


Spin-down luminosity (xlO 38 erg/s) 


4.6 


1.5 



(Finley et al, 1993), GINGA (Deeter et al., 1999), 
poSAX (Mineo et al., 1999), ASCA (Hirayama et al., 2002), 
RXTE (de Plaa et al., 2003), and Chandra (Gotthelf & 
Wang, 2000; Hwang et al., 2001; Kaaret et al., 2001; Petre et 
al., 2007). This source is not detected in 7-rays (> 50 MeV), 
with only upper limits reported by EGRET (Thompson et 
al., 1994). 

In the optical/near UV band PSR B0540-69 is a faint 
~22 magnitude source with a power law spectrum having a 
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steep spectral index a v ~ 1 (Serafimovich et al., 2004) while 
the Crab has a flat spectrum, a v ~ 0, in the same band. 

In the X-rays, this pulsar has a flux of a few milliCrab, 
and a spectrum fitted by various authors as a power law 
with photon indices in the range from 1.8 to 2.0. De Plaa 
et al. (2003), using high resolution RXTE spectra, showed 
evidence for a curved, log-parabolic spectrum, i.e. with the 
form F(E) = KE~ a ' bhosE . A similar spectrum has been 
observed also for the Crab pulsar (Massaro et al., 2000, 2006) 
and few other sources (Cusumano et al., 2001; Mineo et al., 
2004). 

PSR B0540-69's PWN (also called N158A) has a torus- 
like structure, with a major radius of ~4" and with the 
presence of jets along the direction of proper motion. At 
variance with the Crab, the PWN of PSR B0540-69 is fur- 
ther surrounded by a Supernova Remnant (SNR) thermal 
shell, with a radius of about 30", that has radio and soft X- 
ray emission. Optical (HST; Serafimovich et al., 2004) and 
X-ray (Chandra; Hwang et al., 2001; Kaaret et al., 2001; 
Petre et al., 2007) observations showed that the various re- 
gions of the PWN/SNR system have different brightnesses 
and spectra, not described by simple power law models. 

A first analysis of PSR B0540-69 using INTEGRAL 
data was performed by Slowikowska et al. (2007), who re- 
ported pulse profiles in different energy bands and a total, 
phase-averaged, spatially-extended spectrum. 

In this paper we present a comprehensive X/7-ray study 
of this pulsar-PWN system, analyzing in detail the INTE- 
GRAL observations together with new SWIFT/XRT data. 
We also exploit the very large RXTE observation dataset of 
this pulsar, which spans a period of more than ten years. 

This paper is structured as follows. In Sect. [2] we de- 
scribe the X-ray observations and data reduction and we 
present our results on the pulsar and PWN spectra. In Sect. 
[3] and [4] we discuss these results and their interpretation also 
in the context of pulsar-PWN high-energy emission models. 



2 DATA REDUCTION AND SPECTRAL 
ANALYSIS 

2.1 Swift XRT 

PSR B0540-69 was observed fifteen times in three different 
periods by the XRT (X-Ray Telescope) instrument onboard 
the Swift satellite (Gehrels et al. 2004). In Table[2]we report 
the log of these observations, performed between 2005 and 
2007. XRT is a focusing telescope with a CCD detector, and 
its useful band of sensitivity is roughly 0.7-7 keV. The main 
observing modes are Photon Counting (PC) and Windowed 
Timing (WT). In PC mode the data are integrated in a 2.5 
seconds exposure, retaining full imaging and spectroscopic 
capabilities but with very limited time resolution; while in 
WT mode the central 200 CCD columns (8' of the field of 
view) are continuously registered in one-dimensional mode, 
with a time resolution of 1.7 ms. 

Raw event data (level 1) were reduced following 
the standard pipeline (xrtpipeline v. 0.11. 6 task) using 
the HEADAS v.6.4 package, including SWIFT software 
v.2.8 and the latest CALDB calibration files. Events were 
screened adopting standard grade filtering. The total obser- 
vation time is respectively 90.4 ks for the PC mode and 23.5 



Table 2. Log of XRT observations of PSR B0540-69. 



No. 


ObsID 


Start date 


Mode 


Exposure (s) 


01 


00053400001 


2005 Apr 14 


PC 


8 415 


02 


00053400002 


2005 Apr 19 


PC 


3 319 


03 


00053400003 


2005 Apr 20 


PC 


6 644 


04 


00053400004 


2005 Apr 28 


PC 


7340 


05 


00053400005 


2005 May 04 


PC 


8 071 


06 


00053402001 


2006 Jan 26 


WT 


2 815 


07 


00053402002 


2006 Feb 28 


WT 


1505 


08 


00053402003 


2006 Mar 04 


WT 


16 548 


09a 


00053402004 


2006 Apr 18 


WT 


2 651 


09b 






PC 


3 866 


10 


00053402005 


2006 Apr 19 


PC 


11205 


11 


00053402006 


2006 Apr 20 


PC 


6 261 


12 


00053402007 


2007 Nov 05 


PC 


5147 


13 


00053402008 


2007 Nov 08 


PC 


3 260 


14 


00053402009 


2007 Nov 09 


PC 


11431 


15 


00053402010 


2007 Nov 11 


PC 


15 438 



ks for WT mode. We used the PC mode observations to ob- 
tain a total (pulsar+PWN) source spectrum, and the WT 
mode observation to derive the pulsed spectrum. 

2.1.1 Photon Counting mode observations 

Source spectra were extracted from a circular region of ra- 
dius 20 pixels (equal to 47". 2) excluding the central pixel 
to avoid count rates higher than 0.6 cts/s and thus pile-up 
contamination. The PWN radius (~4") is comparable to the 
XRT point spread function, so this source is almost point- 
like. Background spectra were extracted from an annular 
region outside the source region, with radii of 20 and 50 pix- 
els respectively. The corresponding Ancillary Response Files 
(ARF) were generated by means of the xrtmkarf task, tak- 
ing into account the exposure map generated by the pipeline. 

All spectra, rebinned to have at least 20 counts for 
each energy channel, were fitted individually with XSPEC v. 
11.3.1, using an absorbed powerlaw (wabs*powerlaw) model. 
The results were verified to be consistent, and all were found 
to match within the errors: the spectra were then summed 
together. The result is shown in Fig.[T] and the best fit model 
has an hydrogen column density Nh = (3.7±0.1)T0 21 cm -2 
and a photon index T = 1.98 ± 0.02. The reduced x 2 is 1.01 
for 483 degrees of freedom. The 2-10 keV absorbed flux is 
(2.61 ± 0.03) • 10~ n erg cm -2 s" 1 . 

It should be noted that our fit assumes Galactic ele- 
mental abundancies and doesn't discriminate between Milky 
Way and LMC absorption, although the former is about 
one order of magnitude less than the latter. The derived 
Nh value, therefore, should be regarded as a simple, phe- 
nomenological parameter for the description of the spec- 
trum. See Serafimovich et al. (2004) for an extensive discus- 
sion on the photoelectric absorption in the X-ray spectrum 
of PSR B0540-69. 

Our results are in very good agreement with previous 
observations: BeppoSAX LECS and MECS (Mineo et al. 
1999) data gave in fact N H = (3.7±0.1) ■ 10 21 cm -2 and T = 
1.94 ± 0.03, also in agreement with previous ROSAT results 
(Finley et al., 1993). Moreover, summing the Chandra pulsar 
phase-averaged and PWN spectra determined by Kaaret et 
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Table 3. Swift-XRT best fit folding frequencies. 

Obs. # uo (Hz) Reference epoch (MJD) 



channel energy (keV) 

Figure 1. Swift XRT X-ray (0.7-7.14 keV) spectrum of the total 
(pulsar+PWN) emission from PSR B0540-69. 



al. (2001), we also obtain a very good agreement with our 
XRT total spectrum (see Figure [TO]) . 

We also evaluated what is the contribution of the SNR 
thermal shell observed by Chandra on the XRT spectrum. 
From Hwang et al. (2001) it is apparent that, at 20" -40" 
from the pulsar, the SNR falls within XRT spectral extrac- 
tion radius. However, the thermal shell is very soft, with 
almost all of the emission below 2 keV. Moreover, from the 
analysis of Petre et al. (2007), the flux from the shell is only 
a few percent of the total PWN emission. Therefore, the 
SNR contamination on the XRT spectral fit is negligible. 



2.1.2 Windowed Timing mode observations 

In order to obtain the pulse profile and phase resolved spec- 
tra, the photon arrival times of the four WT-mode obser- 
vations (obs. # 06, 07, 08 and 09a, see Table [2j, selecting 
the events in a 40 pixel-wide box centered on the pulsar po- 
sition, were reduced to the solar system barycentre, using 
JE200 ephemeris and the spacecraft orbital data. 

During the observation # 08 the source was erratically 
inside and outside the field of view, therefore considerably 
reducing the useful exposure: we selected only the periods 
with a count rate higher than 0.5 cts/s, obtaining a total 
net exposure time shorter than 2 ks. 

Radio timing observations spanning the XRT data 
epoch are not available, therefore we performed a period 
search, folding the light curves at various trial frequencies 
and computing the corresponding \ 2 value with respect to a 
constant level. Significant pulsation was found only for ob- 
servations # 05 and 09a, the longest ones. The frequency 
corresponding to the maximum \ 2 value was thus obtained 
(Table [3} , and the folded lightcurves were arbitrarily shifted 
to have the pulse maximum (determined by a Gaussian fit) 
at the phase 0.5. The curves were then added together. The 
resulting total light curve, for a net exposure time of 5472 
seconds, is shown in Figure [2] 

A phase- averaged spectrum was computed and yielded 
results that were perfectly consistent with the PC-mode 
observations. We assumed the phases 0.90-1.10 as the off- 
pulse, representative of the nebular spectrum, while 0.10- 
0.90 as the pulsar phase. We subtracted the off-pulse from 



06 19.7609226 
09a 19.7595952 



53761.762 
53843.561 



PSR B0540-69 - Swift/XRT 



Figure 2. Folded lightcurve in 20 phase bins for XRT observa- 
tions of PSR B0540-69, in the energy range from 0.7 to 7.14 keV. 



the pulsed data, obtaining a background-subtracted pulsed 
spectrum, that was fitted with an absorbed power-law, us- 
ing the Nh value derived from the total emission analysis, 
obtaining V — 1.6 ± 0.2 and a 2-10 keV absorbed flux of 
(6.4 ± 0.8) • 10~ 12 erg cm -2 s _1 , in agreement with the pre- 
vious BeppoSAX observations (Mineo et al., 1999). 



2.2 INTEGRAL IBIS/ISGRI 

The INTEGRAL (International Gamma Ray Astrophysics 
Laboratory) mission observed the Large Magellanic Cloud 
during nine orbits in two periods, January 2003 and Jan- 
uary 2004 (Gotz et al. 2006). This satellite has two main in- 
struments, the imager IBIS and the spectrometer SPI. IBIS 
(Imager on Board Integral Satellite, Ubertini et al. 2003) is 
a coded-mask imager and is constitued by two simultane- 
ously operating detector layers, ISGRI and PICsIT. ISGRI 
(Integral Soft Gamma-Ray Imager, Lebrun et al. 2003) is 
an array of 8 CdTe solid state detection units, each having 
32x64 pixels, thus making an overall 128x128 array, that 
has a 12' (FWHM) image resolution and a nominal operat- 
ing range from 15 keV to 1 MeV. Each INTEGRAL orbit 
lasts about three days and consists of several pointings, be- 
cause of the "dithering" observational mode. Each pointing, 
the duration of which is typically half an hour, forms a "Sci- 
ence Window" (ScW) . 

PSR B0540-69 was observed during revolutions 27, 28, 
29, 33, 34, 35 (January 2003) and 150, 151, 152 (January 
2004). Individual pointings for all the observations per- 
formed by the first 3.5 years of INTEGRAL observations 
were processed using the standard software package OSA 
v.5.1 (Goldwurm et al., 2003) and mosaicked to obtain an 
all-sky image. PSR B0540-69 is then detected at 6.5<r sig- 
nificance for about 221 ks of corrected on-source exposure 



4 R. Campana et al. 



PSR B0540-69 - INTEGRAL/IBIS-ISGRI 



channel energy (keV) 

Figure 3. INTEGRAL IBIS-ISGRI (20-100 keV) spectrum of 
PSR B0540-69. 
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Figure 5. Eleven years of PSR B0540-69 RXTE-PCA timing. In 
the lower panel are reported the residuals after the fit with Eq. 
0, which is shown as the red continuous line in the upper panel. 



PSR B0540-69 - Swift/XRT + INTEGRAL/IBIS-ISGRI 




1 2 5 10 20 50 



channel energy (keV) 

Figure 4. Swift XRT (0.7-7.14 keV) and INTEGRAL IBIS- 
ISGRI (20-100 keV) spectrum of PSR B0540-69. 



(see Bird et al., 2007, for the details of data analysis). The 
spectrum is obtained calculating the mean flux in various 
energy bands spanning the range 20-100 keV. 

A power-law fit of IBIS data in the 20-100 keV band is 
reported in Fig. [3] The photon index is T = 2.12 ± 0.25 and 
the reduced X 2 is 1-5 (for 5 d.o.f.). The 20-100 keV flux is 
2.9 jl^'j ' 10 -11 erg cm -2 s - , in agreement with Gotz et al. 
(2006). 

Note that this photon index is steeper than the XRT 
one, although compatible because of the large uncertainties. 
In fact, a simultaneous fit of XRT and IBIS data for the 
total source (PSR+PWN) spectrum is reported in Fig. [4] 
The resulting parameters are: Nh = (3.8 ±0.1) ■ 10 21 cm -2 , 
F = 1.97 ± 0.02 and X 2 = 1-09/489 d.o.f. The intercalibra- 
tion factor between IBIS and XRT data is / = 1.05 ± 0.08, 
and clearly the fit is constrained by the high-quality XRT 
data. The IBIS spectrum thus can be regarded as a smooth 
extrapolation of the XRT data. 



PSR B0540-69 - RXTE/PCA 




0.5 1 1.5 2 

Phase 



Figure 6. Total light curve obtained from the full PCA dataset, 
in the energy range 2—60 keV, with 100 phase bins. 

2.3 RXTE 

2.3.1 PCA 

PSR B0540-69 was observed several times with the Rossi 
X-Ray Timing Explorer (RXTE) satellite. We used data ob- 
tained with the Proportional Counter Array (PCA, Jahoda 
et al., 1996) instrument in the Good Xenon mode, that en- 
sures the full timing accuracy (about 1 ^ts). PCA is a 1° x 1° 
collimated array of proportional counters, and is composed 
of 5 PCUs (Proportional Counter Units). The useful sensi- 
tivity range is about 2-60 keV. 

For the timing analysis we selected all the archive ob- 
servations having PSR B0540-69 at less than 30' from the 
centre of the field of view (see Table [4} and with a duration 
longer than 3000 seconds. The data set spans all years time 
range, from 1996 to 2007. The observations were screened 
for South Atlantic Anomaly passages, elevation angle and 
electron background, following the standard criteria, and 
barycentred using the f axbary task of the FTOOLS pack- 
age. The total screened exposure time is more than 2 Ms, 
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Table 4. Log of RXTE observations of PSR B0540-69. Each proposal ID contains several observations (from 2 to 95). The offset is 
between the center of the pointing field of view and the location of the pulsar. The on-source exposure time is after the screening (PCA) 
and dead-time (HEXTE) corrections. 



Proposal 




First- 


Last 






Unset 


rl/A 


T_T T7 V T7 A 




1JJ 




Observation Date 




(arcmin) 


exposure (s) 


exposure (s) 


exposure (s) 


10206 


1996 


Aug 11 - 


1996 


Nov 


17 


0.30 


54 049 


35 200 


34 955 


10218 


1996 


Oct 12 - 


1996 


Dec 


22 


25.0 


92 086 


29 306 


29 212 


10250 


1996 


Feb 02 - 


1996 


Oct 


04 


24.6 


37 837 


12 422 


12 325 


20188 


1996 


Nov 30 - 


1997 


Dec 


12 


24.5 


301 275 


118 935 


118114 


30087 


1998 


Jan 04 - 


1998 


Sep 


30 


24.6 


103 455 


44 566 


44 265 


40139 


1999 


Jan 19 - 


2000 


Feb 


15 


15.8 


185 520 


71834 


71 439 


50103 


2000 


Mar 13 - 


2001 


Mar 


15 


15.8 


219 360 


178 802 


177 325 


50414 


2000 


Jun 05 - 


2000 


Jun 


23 


0.30 


5 552 


3 634 


3 601 


60082 


2001 


Mar 26 - 


- 2003 


Feb 


21 


15.8 


281719 


223 510 


221 934 


70092 


2002 


Mar 21 - 


2003 


Mar 


30 


15.8 


268 534 


233 883 


232 737 


80089 


2003 


Apr 16 - 


2004 


Oct 


10 


15.8 


263 151 


211649 


210143 


80118 


2004 


Jan 07 - 


2004 


Jan 


12 


24.6 


40 977 


21593 


21606 


90075 


2004 


Mar 03 - 


2005 


Mar 


29 


15.8 


83 507 


67 633 


75 325 


91060 


2005 


Jun 20 - 


2006 


May 


28 


15.8 


47523 


32 707 


39 855 


92010 


2006 


Apr 25 - 


2007 


Jun 


22 


15.8 


33145 


27210 


29 755 


Total (Ms) 














2.018 


1.313 


1.322 



a significant improvement, by a factor of ~3, with respect 
to the previous RXTE analysis of de Plaa et al. (2003), 
which had a total exposure of 684 ks. A period search was 
performed for each observation, and the found frequencies 
(Figure O were fitted with a second-order Taylor expansion 
of the frequency history, 

v(t) = v + i>o(t - to) + y(£-M 2 (1) 

in order to find a satisfactory ephemeris spanning the data 
set. The values derived from the fit are, for the time range 
in MJD 50100-54100: u = 19.774 Hz, u = -1.87253- 10 -10 
Hz/s and vq = 3.69 ■ 10~ 21 Hz/s 2 , for the reference epoch to 
MJD 52954. In the 11-years time interval of RXTE observa- 
tions the spin frequency of the pulsar decreases from a value 
of about 19.82 Hz to 19.75 Hz. Note that the XRT folding 
frequencies (Sect. [2~1.2|l are in good agreement with the ex- 
trapolation of these ephemeris at the relative epochs. The 
values of the frequency and its time derivative are in good 
agreement with other works (Cusumano et al., 2003; Living- 
stone et al., 2005; Slowikowska et al., 2007). However, the 
value of the second time derivative is not well constrained by 
the fit. Each observation was thus phase-folded. The align- 
ment between the pulse profiles of each observation was per- 
formed by a cross-correlation method, first with a sinusoidal 
function in order to have the peak maximum at phase 0.5, 
then with an analytic template based on the resulting total 
light curve. A 100-bin pulse profile in the energy band 2-60 
keV is showed in Figure [6] and energy-resolved profiles are 
showed in Figure [7] 

From this figure it is apparent that there is no evident 
change of the light curve shape with energy. Following de 
Plaa et al. (2003), we fitted the full pulse profile (Figure 
[6]) with an analytical model described by the sum of two 
Gaussians and a constant. We then performed this fit on 
all the Figure [7] light curves, fixing the central phase and 
the widths of the Gaussians at the values determined in the 
total light curve fit, and we found no statistical significant 
variation on the relative normalizations of the Gaussians. We 



also checked for the stability of the pulse profile, comparing 
it as observed at various epochs, and we did not found any 
significant changes in shape or flux. 

For the spectral extraction we used only the observa- 
tions pertaining to the last gain epoch (Epoch 5, from May 
13 2000 onward, proposal ID 50103 to 92010) in order to 
ensure the consistency of the analysis. The resulting dataset 
has a total exposure of 1.2 Ms. In a similar way with respect 
to XRT, the phases 0.10-0.90 were selected as the pulsed sig- 
nal, and the phases 0.90-1.10 as the offpulse, that was sub- 
tracted on the pulsed signal. Offset-corrected response ma- 
trices were determined for each observation and each PCU, 
and added together weighting them by their exposure. Spec- 
tra were also extracted for each PCU, rebinned to have a 
minimum of 20 counts for channel (however, the rebinning 
affects only the high energy channels, over 30 keV), and ver- 
ified for consistency. 

A power-law fit of the pulsed signal, in the energy range 
2-30 keV and with Nh frozen at the value given by the XRT 
analysis, gave the values F = 1.89 ± 0.01 and % 2 = 1-14/59 
d.o.f.; the residuals, however, showed a systematic behaviour 
typical of curved spectra, with an excess in the central bins 
(see Figure |8j left panel). We considered therefore a log- 
parabolic law, and the fit gave the values a = 1.47 ± 0.01, 
b = 0.25 ± 0.06 and \ 2 = 0.96/58 d.o.f., and the residuals 
don't show any systematic deviation (Figure [8] right panel). 
We performed an F-test in order to assess the statistical 
significance of the log-parabolic model versus the power law, 
and the null hypothesis probability resulted 9.7 ■ 10~ 4 . The 
2-10 keV pulsed absorbed flux is 6.5± ;7-10~ 12 erg cm" 2 s" 1 , 
in agreement with the XRT value. 

The 30-60 keV PCA spectrum has a much lower statis- 
tics. Extending the previous fit up to 60 keV does not alter 
the values of the parameters. The 30-60 keV spectrum alone 
is well described by a single power law with a photon index 
around 2. 

There are some discrepancies with respect to the spec- 
tral results of de Plaa et al. (2003), both in the curvature 
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i i i i i i i i i i i i — i i i i i i i i i i i — i i i i i i i i i i i i 

0.5 1 1.5 0.5 1 1.5 0.5 1 1.5 2 
Phase Phase Phase 

Figure 7. Light curves in various PCA energy bands, from 2 to 60 keV. The curves are with 50 phase bins, except for the last two that 
have 25 bins. 




parameter and in flux. Firstly, they gave the log-parabolic 
law in terms of natural, instead of decimal, logarithm. How- 
ever, also correcting for the change of logarithm basis, their 
curvature parameter is larger (i.e. the log-parabola is more 
curved) than ours. There are many possible explanations 
for the discrepancy: (a) the different methods used to deter- 
minate the pulsed spectrum (de Plaa et al. determined the 
pulsed counts excess in 14 energy bins by means of an ana- 



lytical fit of the light curve); (b) the different energy range 
of their fit (that included also soft X-ray ROSAT data), 
together with a different Nh value, could affect the curva- 
ture; (c) their spectral PCA dataset has an exposure (462 
ks) much smaller than ours (1.2 Ms); and (d) a real change 
in the spectrum between the calibration epoch 3 (1996-99) 
dataset used by de Plaa et al. and our epoch 5 dataset (2000- 
07). Regarding the first of the above points, we determined 
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the excess counts in various energy bands by means of ana- 
lytical fits and we found them in agreement with the pulsed 
counts determined by our offpulse-subtraction method. We 
also performed a spectral analysis of the epoch 3 data (pro- 
posal ID 10206 to 30087) with our approach, and we found 
that both the spectral parameters and the flux are perfectly 
la compatible with the better statistics of the epoch 5 spec- 
trum. Moreover, the flux values of de Plaa et al. are greater 
than ours by a factor of about 30%. This discrepancy, too 
great to be due only to the different estimation of the pulsed 
counts above the offpulse level, could be due to the different 
way of estimating the effective exposure time. We verified, 
in fact, that our data agree very well with the ones of de 
Plaa et al. when the exposure is reduced by a factor equal 
to the pulsed phase fraction. Note that the Chandra pulsed 
spectrum in the range 0.6-10 keV (Kaaret et al., 2001, also 
reported by Serafimovich et al., 2004) matches our data, as 
seen in Figure [TU] 



2.3.2 HEXTE 

Data from the instrument HEXTE (High Energy Timing 
Experiment, Rothschild et al., 1998) onboard RXTE were 
also used. HEXTE employs NaI(Tl)-CsI(Na) phoswich scin- 
tillation detectors, and consists of two separated clusters 
(A/0 and B/l), each containing four detectors. Detector 2 
of cluster B has lost its spectral capabilities, and is useful 
only for timing analyses. We extracted the data using the 
science event mode that ensures the full timing resolution 
(8 lis), and we barycentered and folded the data using the 
PCA ephemeris. The total on-source, dead-time corrected 
exposure is about 1.3 Ms for each cluster. 

Phase alignment was performed using the shift deter- 
mined from the simultaneous PCA observations. The folded 
lightcurve is showed in Figure [9] To increase the signal to 
noise ratio for the spectral extraction, we used only the ob- 
servations having an offset less than 16', being the response 
dramatically dependent on the offset angle. We produced the 
offset-corrected response matrices, using the task hxtrsp, 
and added them together weighting them by the exposure, 
in a similar way to what was done for PCA. The pulsed 
spectrum was extracted in the same phase bins of PCA. In 
the 15-250 keV range it is well approximated by a single 
power law with a photon index of 2.1 ± 0.2. The 20-100 
keV pulsed flux is &.lt° 2 \ ■ 10~ 12 erg cm" 2 s" 1 . Note that 
this photon index, although affected by a rather large uncer- 
tainty, is higher than the one measured in the keV range, and 
shows a high energy spectral steepening. We see (Fig. llQ[) , in 
fact, that the HEXTE spectrum is fully compatible with the 
log-parabolic fit of PCA data, and follows smoothly the cur- 
vature up to 100 keV. HEXTE data, on the contrary, would 
lie systematically below the extrapolation of the power law 
fit of PCA data. 



3 THE X-RAY SPECTRA OF THE PULSAR 
AND THE PWN 

The spectral analysis of the total (pulsar + PWN) and 
pulsed emission from PSR B0540-69, described in the previ- 
ous Section, is summarized in Table [5] The Spectral Energy 
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Figure 9. Total light curve obtained from the HEXTE dataset, 
in the energy range 15-100 keV, with 20 phase bins. 



Distributions of the total and pulsed emission are plotted in 
Figure [T0l 

The agreement between the various instruments in the 
1-10 keV band is very good, both for the total emission and 
the pulsed one. 

It's clearly visible how the pulsed emission accounts 
only for about 20% of the total emission in the soft and 
hard X-rays. The ratio of the total source flux to the pulsed 
flux is in fact ~ 5 in the 20-100 keV band (ISGRI total 
versus HEXTE pulsed spectrum) and ~ 4 in the 2-10 keV 
band (XRT total versus PCA pulsed spectrum). 

Therefore, subtracting the pulsed flux from the to- 
tal one to evaluate the PWN emission in the range 1-100 
keV, we obtain an energy- dependent flux ratio Jpwn /Fpsr 
between 2.5 around the peak of the pulsed emission at 
about 12 keV, and ~6 at the range boundary. Integrating 
over the whole energy range, we derive a luminosity ratio 
£pwn/£psr ~ 3.6. This value is expected to correspond 
also to the ratio between the efficiencies 7?pwn/?7psr of con- 
version of the spin-down luminosity E to electromagnetic 
radiation, and is in agreement with the mean value (~4) for 
a large sample of X-ray pulsars (Kargaltsev & Pavlov, 2008) 

For the total energy conversion efficiency, using E = 
1.5 • 10 38 erg/s, we obtain a value of tjtot = 5.8% for the 
20-100 keV energy range. The PWN emission contributes 
to this figure for the 4.6%, while the pulsed emission for the 
1.2%. Moreover, we have t/tot = 5.2% for the 2-10 keV 
energy range (PWN: 3.9%, pulsar 1.3%). 

If we subtract the log-parabolic, pulsed spectrum to the 
power-law, total spectrum, we can derive an approximate 
photon index for the nebular spectrum. For the 2-10 keV 
band we obtain F ~ 2.05, slightly softer than that of the 
total emission, in agreement with the mean photon index of 
Chandra nebular and offpulse spectra (~ 2.04, Kaaret et al., 
2001). 

Petre et al. (2007) have shown that in the 1-10 keV 
band the PWN photon index increases from a value of ~1.4 
near the pulsar, to ~2.5 at a radius of 4". The nebular size, 
however, shrinks with the energy. An extrapolation of the 
relation found by Petre et al. to the IBIS energy range re- 
sults in a PWN radius of about 2.5": inside this radius, the 
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Table 5. Results of the spectral analysis. The absorbed flux values are in units of fO 11 erg cm 2 s 1 . 



Instrument 


Parameter 


Flux (2-10 keV) 


Flux (20-100 keV) 


Swift/XRT - Total 


r = 


1.98 ±0.02 


2.61 ± 0.03 




INTEGRAL/IBIS - Total 


r = 


2.12 ± 0.25 




2 9+ ai 


Swift/XRT - Pulsed 


r 


= 1.6 ± 0.2 


0.64 ±0.08 




RXTE/PCA - Pulsed 


o = 

6 = 


1.47 ±0.01 
0.25 ±0.06 


65+ 04 
U - DO -0.07 




RXTE/HEXTE - Pulsed 


r 


= 2.1 ± 0.2 




61+ - 06 



" ■ INTEGRAL IBIS/ISGRI - Total 

- • SWIFT XRT - Total 

- * SWIFT XRT - Nebula 
♦ SWIFT XRT - Pulsed 
■ RXTE PCA - Pulsed 

- • RXTE HEXTE - Pulsed 




Energy (keV) 



Figure 10. Spectral Energy Distribution (SED) for the emission from PSR B0540-69 and its PWN. The red continuous line is the 
log-parabolic law that fits PCA pulsed data (see Sect. 12.3.11 . while the black dotted line is the power-law simultaneous fit of XRT and 
IBIS total data (see Sect. [2T2t . The black dotted and red dashed areas are the ±1<t Chandra pulsed and total (pulsar phase averaged + 
PWN) data, respectively (from Kaaret et al., 2001, and Serafimovich et al., 2004). 



mean photon index is about 1.8. Again, this is consistent 
with the derived PWN spectrum obtained subtracting the 
pulsed contribution to the total spectrum. For the 20-100 
keV band, being after the peak of the pulsed spectrum, the 
PWN has an harder spectrum than the total emission, with 
a photon index around 1.8. This qualitative agreement, how- 
ever, relies on the assumptions that both the nebular photon 
indices and the parameterization of the nebular radius could 
be extrapolated to the hard X-rays, and of course it is af- 
fected by the low statistics data for both the total and the 
pulsed emission. 

It is interesting to compare the soft 7-ray properties 
of PSR B0540-69 with the other INTEGRAL-detected pul- 
sar/PWN systems (PSR J1846-0258, McBride et al., 2008; 
PSR J1617-5055, Landi et al., 2007; Vela, Hoffmann et al., 
2006; PSR J1513-5906, Forot et al., 2006; PSR J1811-1925, 
Dean et al., 2008). With r\ ~ 5.8% for the total emission, 
this source has the highest conversion efficiency: the mean 
of the sample is around 1%, while Vela is a notable excep- 



tion with rj ~ 0.02%. Moreover, the ~2.1 photon index of 
PSR B0540-69 fits squarely in the observed range for this 
sample of IBIS sources (2.0-2.3). The majority of these sys- 
tems have the hard X-ray emission dominated by the PWN, 
except for PSR J1617-5055 and PSR J1513-5906. The rel- 
ative contribution from the pulsar and the PWN, however, 
differs very much from source to source. A more detailed 
comparison will be the subject of a future paper (Dean et 
al., in preparation). 

Most of the INTEGRAL Galactic PWNs are likely as- 
sociated with HESS TeV sources (e.g. PSR J1617-5055 and 
PSR J1811-1925). However, estabilishing the PSR B0540- 
69 possible nebular TeV emission is difficult, because of the 
low declination of this source and, mainly, of the high dis- 
tance. If we assume that the synchrotron nebula powered 
by PSR B0540-69 has the same intrinsic TeV brightness of 
the Crab, the ground flux is reduced to a value of the order 
of some milliCrab (by a factor ~400). Very long (>100 h) 
dedicated HESS pointings, therefore, are required to detect 
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Figure 11. Distance-corrected luminosity for PSR B0540-69 and Crab. The distance values used are d = 50 kpc and d = 2 kpc, 
respectively. Crab data is from the four NFI onboard BeppoSAX and were taken from Massaro et al. (2006). Also plotted are the 
log-parabolic model for PSR B0540-69 (blue dashed line) and the total multicomponent model for Crab (black continuous line, with the 
Co and Cx components as orange dash-dotted line and violet dash-dot-dotted line). See text for discussion. 



the possible emission (Hofmann, 2001). Another intriguing 
possibility is that this system lacks an appreciable intrinsic 
TeV emission. It is thought (de Jager & Djannati-Ata'i, 2008) 
that efficient high energy 7-ray emission is suppressed when 
the synchrotron losses dominates over the Inverse Compton 
losses, and this can happen, for example, for high magnetic 
energy densites or low background radiation fields, although 
there is no compelling evidence for both hypotheses. 

Petre et al. (2007), under reasonable assumptions, de- 
rive an equipartition magnetic field for the PWN, and found 
a value of B ~ 8 ■ 10~ 4 G, higher than that estimated for the 
Crab by a factor ~2-3 (de Jager & Harding, 1992), and im- 
plying a synchrotron lifetime of electrons of about 3 years. 
The efficiency for the very high energy IC emission should 
therefore be lower with respect to the Crab, making further 
difficult the detectability of TeV flux. It will be interesting 
to see if GLAST will be able to detect nebular emission in 
the GeV range, and to check if the derived PWN parameters 
are consistent with the previous estimates. 

Note also that an EGRET source, 3EG J5033-6916, 
is located about 0.6° south-east from the pulsar (thus 
compatible with the PSF radius), but its association with 
PSR B0540-69 is unlikely. This 7-ray source has been in- 
terpreted as due to diffuse emission originating from the 
interaction between cosmic rays and the LMC interstellar 
medium (Sreekumar et al., 1992; Lin et al., 1996), and the 
observed flux is consistent with the theoretical predictions 
(Fichtel et al., 1991). 



3.1 Pulsed spectrum 

Massaro et al. (2006) developed a phenomenological model 
of the Crab pulsar X and 7-ray emission, explained as due to 
two different emission components Co and Cx , having a log- 
parabolic spectrum with the same curvature (b ~ 0.16) but 
different energy peaks and phase distributions. This multi- 
component (MC) model also describes the energy variations 
of the pulse profile. The 7-ray emission above 10 MeV is 
well reproduced, both in the total and phase-resolved spec- 
tra, assuming that the Co and Cx have an higher-energy 
counterpart Co 7 and Cx-y, with the same phase distribu- 
tions and the same log-parabolic spectra at different peak 
energies. 

In Figure ITT1 the pulsed spectra of both PSR B0540-69 
and Crab are plotted, rescaled for the distance values (we 
assumed d = 50 kpc and d — 2 kpc, respectively), together 
with the Co and Cx components of the Crab MC model. 
The emission of PSR B0540-69 is higher by a factor of about 
2 than the Crab, although the true value is strongly depen- 
dent on the assumed distance (see Schaefer et al., 2008, for a 
discussion on the LMC distance value), but also on the dif- 
ference in the solid angle subtended by the emission beam, 
here assumed to be isotropic. If we adopt the distance values 
from HSTKP (50.1 ± 2.4 kpc, Freedman et al., 2001), the 
luminosity ratio lies in the range of 1.1-3.7. 

It is interesting to note that, although PSR B0540-69 
spectral curvature is larger (b ~ 0.24 versus b ~ 0.16 for the 
Crab), the peak energy is sensibly close to the E p ~ 12 keV 
value of the Co component of the MC model. The absence 
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of a variation with energy of the light curve also makes plau- 
sible the interpretation that the PSR B0540-69 emission is 
only due to its counterpart of the Crab Co component. The 
pulse shapes of these pulsars are however very different, and 
this could be due either to different viewing conditions or to 
unlike properties of the emission regions. 

Takata & Chang (2007) recently developed an outer- 
gap model of Crab and PSR B0540-69 emission, extending 
the emission region also inside the null-charge surface, along 
the last closed field line. Their results are in broad agree- 
ment with the MC phenomenological model, if we identify 
the Co component as due to the emission beyond the null 
charge surface to the light cylinder, and the Cx component 
as coming from the region below the null charge surface 
up to the neutron star surface. In their model, the X-ray 
spectrum is due to synchrotron emission from secondary e 
pairs, and has a curved shape similar to our log-parabolic fit. 
Takata & Chang (2007) explain the different pulse shapes 
of Crab and PSR B0540-69 assuming a different inclination 
angle (a ~ 50° and a ~ 30°, respectively) between the spin 
axis and the magnetic axis, and a different viewing angle £. 
Also, the latter pulsar has a thicker emission region of the 
secondary pairs. 

If we assume that also PSR B0540-69 has a Co-, com- 
ponent at the same peak energy (300 MeV) of the Crab 
equivalent, and with the same relative normalizations of the 
two log-parabolic spectra, we can obtain a value of the flux 
at the peak frequency: ~ 4 • 10~ 7 MeV cm -2 s _1 . This 7- 
ray flux prediction is of the same order of magnitude with 
respect to the one derived from Takata & Chang (2007) by 
Inverse Compton scattering. The predicted pulsed flux is 
about one order of magnitude lower than the 5a one-year 
GLAST survey sensitivitjfl Note that the EGRET upper 
limit (Thompson et al., 1994) is about I0~ 5 MeV cm" 2 s" 1 
in the energy range above 100 MeV. 



4 CONCLUSIONS 

PSR B0540-69 is an interesting young pulsar, showing many 
similarities but also some differences with respect to the 
Crab. The study of this pulsar, although made difficult by its 
high distance and apparent faintness, could help the theoret- 
ical understanding of the characteristics of the high-energy 
emission from the pulsars and their PWNs. 

With the analysis of the complete RXTE dataset of 
observations of PSR B0540-69, exploiting for the spectral 
analysis the longest gain epoch that spans more than half 
of these data, together with that of new SWIFT/XRT and 
INTEGRAL/IBIS data, we have compiled an accurate and 
up-to-date picture of the emission from this source in the 
X and soft 7-rays, from 0.7 to 200 keV. Thanks to the full 
HEXTE dataset we have also determined the contribution 
of pulsed hard X-ray emission to the total one. 

Regarding the pulsed emission, we confirm the pres- 
ence of a curved spectrum, peaking at about 12 keV, well 
described by a log-parabolic law up to >100 keV. We have 
shown also that it presents some similarities with the Crab 
pulsed emission. 



1 http://www-gl.ast . slac . Stanford. edu/sof tware/IS/glast_lat 



We have studied in detail the total (pulsar+PWN) 
emission from this source, in the energy range 0.7-100 keV, 
and we have shown that the major (about 75-80%) con- 
tribution comes from the synchrotron nebula. The compar- 
ison of PSR B0540-69 with other IBIS pulsar/PWN sys- 
tems shows that it has the highest conversion efficiency in 
the hard X/soft 7-ray band, with similar spectral proper- 
ties with respect to the other sources. Most of these INTE- 
GRAL PSR/PWN systems have been associated with HESS 
sources. PSR B0540-69 looks different, in this respect, not 
showing evidence of emission in TeV range. We have dis- 
cussed the possible instrumental limitations, due to the low 
flux and declination, but there is also the possibility that 
this source genuinely lacks TeV emission. 

Future observations, especially in the hard X-ray band, 
will further refine this picture. High angular resolution long 
Chandra X-ray observations will permit an accurate study 
of the synchrotron nebula to derive its physical parameters. 
Additional INTEGRAL pointings will increase the statistics 
of the total emission, while the advent of hard X-ray focusing 
telescopes, like Simbol-X (Ferrando et al., 2005), will help to 
disentangle the pulsed emission from the nebular one also in 
the HEXTE and IBIS energy range. Moreover, deep GLAST 
observations could be useful to estabilish more strict upper 
limits to the 7-ray emission from this source. 
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